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ABSTRACT 

We present a new photometric catalog of 326 candidate globular clusters (GCs) in the nearby spiral 
galaxy MIOI, selected from B, V, and I Hubble Space Telescope Advanced Camera for Surveys 
images. The luminosity function (LE) of these clusters has an unusually large number of faint sources 
compared with GCLEs in many other spiral galaxies. Accordingly, we separate and compare the 
properties of “bright” {My < —6.5) versus “faint” {My > —6.5; one magnitude fainter than the 
expected GC peak) clusters within our sample. The LF of the bright clusters is well fit by a peaked 
distribution similar to those observed in the Milky Way (MW) and other galaxies. These bright clusters 
also have similar size (rgff) and spatial distributions as MW GCs. The LF of the faint clusters, on 
the other hand, is well described by a power law, dN{Ly)/dLy oc Ly with a = —2.6± 0.3, similar to 
those observed for young and intermediate-age cluster systems in star forming galaxies. We find that 
the faint clusters have larger typical Ceff than the bright clusters, and have a flatter surface density 
profile, being more evenly distributed, as we would expect for clusters associated with the disk. We 
use the shape of the LF and predictions for mass-loss driven by two-body relaxation to constrain the 
ages of the faint clusters. Our results are consistent with two populations of old star clusters in MlOl: 
a bright population of halo clusters and a fainter, possibly younger, population of old disk clusters. 
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1. INTRODUCTION 

Ancient (> 10 Gyr) star clusters formed during the 
early assembly of most galaxies, and therefore give in¬ 
sight into the broad formation history of their hosts. El¬ 
liptical and lenticular galaxies host red, metal-rich glob¬ 
ular clusters (GCs) believed to be associated with their 
bulges, and blue, metal-poor GCs believed to be associ¬ 
ated with their halos. It is not yet clear whether these 
different populations have different ages in all early-type 
gala xies, but an age d ifference has been observed in sev¬ 
eral (|Park et al.ll2012l) . Spiral galaxies also form metal- 
poor halo and metal-rich bulge/thick disk clusters, al¬ 
though the fraction of red-to-blue GCs is typical ly lower 
than f ound in similar m ass early-type galaxies. [Minnitil 
(|1995f) and iCotd (|1999r) showed that metal-rich Milky 
Way (MW) GCs in the inner regions of the Galaxy are 
associated with the bulge. 

In addition to ancient GCs, spiral galaxies form 
younger clusters with a large range of ages in their disks. 
In the MW, young (< 100 Myr) open clusters are con¬ 
fined to the thin disk an d old (> 3 Gyr) open cluster s are 
foun d in the thick disk (iPort egies Zwart et al.ll20I0[ ). In 
fact, iKharchenko et aP (j20I3l l find a smoothly increas¬ 
ing dispersion in the distance of MW star clusters from 
the Galactic plane with increasing cluster age (see their 
Figure 5). 

The GC luminosity functions (LFs) of both elliptical 
and spiral galaxies, including the MW, are observed to 
have a peaked shape with a “turnover” caused by the 
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earlier disruption of lower mass clusters due to “evapo- 
ration” of stars by two-body relaxation in a tidal field 
(iFall fc Zhand[200l[l . Small variations in the peak lumi¬ 
nosity and width of GCLFs are believed to come from 
possible differences in cluster ages, cluster densiti es, or 
the strength of the tidal field (|Villegas et ahllMT^ . De¬ 
spite the known small variations, the turnover of GCLFs 
is nearly universal, which has led to its use as a standard 
candle for distance determinations (|Harrisll200lL p. 223). 

It is therefore important to study cases where the 
GCLF does not follow the “universal” pattern. An in¬ 
teresting example is the so-called “faint fuzzies” that 
have been foun d in the lenticular ga l axies NGC 1023 
and N GC 3384 ((Brodie fc Larsenll2nn2l : iLarsen fc Brodiel 
l2000tl. lenticular and ell iptical galaxies in the Virgo clus- 
ter (IPeng et al.l l2006bll. and M51 and its companion 
(iScheepmaker et al.l 120071: iHwang fc Led 120061) . These 
clusters are quite extended with effective radii (res) of 
7-15 pc and ages upwards of 7-8 Gyr. Despite their in¬ 
termediate/old ages, their LF is observed to continue to 
rise beyond the expected turnover luminosity. Interme¬ 
diate age clusters found in a few early-type galaxies may 
be responsible for discrepancies between the distances 
measured from GCLFs and those mea sured from surface 
brightness fluctuations ()Richtleill2nn3l p. 281). 

Faint fuzzies also have reg more than twice that of typ¬ 
ical G Cs in most ellipicals and spirals. iBrodie fc LarserJ 
(I2002D proposed that perhaps faint fuzzies only form in 
the environment found in lenticuler galaxies, or that per¬ 
haps the y were accreted along with a host dwarf galaxy. 
However. iChie.s-Santos et abl (j2013fl used kinematics and 
spatial comparisons with planetary nebulae and H 1 
to conclude that faint fuzzies may simply be akin to 
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old open disk clusters. They suggest that the reason 
faint fuzzies are observed in lenticulars but not in spi¬ 
rals is that these faint clusters are quite difficult to pick 
out against the strong structure and variable luminosity 
caused by star formation in spiral disks. 

Thus far, most detailed studie s of GCs have fo ¬ 
cused on early type galaxies (e.g. iPeng et ^ (l2006all : 
iVillegas et al.l (1201011 b the MW, M31, and a hand- 
ful of other, m o stly b ulge dominated systems (e.g. 
iGoiidfrooii et al.l (120031 1b Late-type spirals, such as 
MlOl, have no t beeii as thoroughly studied. Previously, 
IChandar et al.l (I2004H studied the GC systems in five spi¬ 
ral galaxies using Hubble Space Telescope (HST) Wide- 
Field Planetary Camera 2 (WFPC2) observations and 
compared with the known distribution of GCs in the 
MW. They found for MlOl and NGC 6946 that the LF 
continues to rise beyond the expected GCLF turnover; 
however, this conclusion was based on a small number 
of cluste rs resulting from the partial coverage of each 
galaxy. iBarmbv et al.l (1200611 studied a larger sample 
of MlOl clusters (1715 clusters with my < 23), but 
they focused on analyzing “blue” clusters {{B — V)o < 
0.45) rather than the older, redder clusters we seek to 
study here. In this paper, we use the same BVI HST 
Advanced Camera f or Surveys (ACS) observations as 
IBarmbv et al.l (|2006ll to examine the LF, colors, sizes, 
and spatial distributions of the red clusters. 

The rest of the paper is organized as follows: we de¬ 
scribe the observations, selection criteria and methods, 
and completeness in §2. In §3, we show the luminos¬ 
ity distribution, color-color plot, color histograms, color 
magnitude diagram (CMD), sizes, and spatial distribu¬ 
tion of our cluster candidates, and we discuss the results 
in §4. Finally in §5, we list our conclusions and describe 
ongoing work on the MlOl cluster system. 

2. OBSERVATIONS, CLUSTER SELECTION, AND 
COMPLETENESS 

2.1. Observations 

Ten pointings within MlOl were taken with the 
HST/ACS Wide Field Channel (WFC) in November 
of 2002 (Program ID: 9490, PL K. Kuntz) using the 
F4351F (B), F555W (F), and F8141F (/) filterfl 
Each field covers 3.3' x 3.4' (see Fig. [T]). Assum¬ 
ing a distance to MlOl of 6.4 ± 0. 2 (ran dom) ±0.5 
(systematic) Mpc (|Shappee fc Stane^ 1201 111 , each field 
covers a 6.1 x 6.3 kpc^ region. Images were pro- 
cess ed through the HLA M ultiDrizzle Pipeline Version 
1.0 (iKoekemoer et al.l 1200311 . which includes bias sub¬ 
traction, cosmic ray rejection by combining two sub¬ 
exposures from CR-Split observing, dark subtraction, 
flat fielding, and drizzling. 

We detect sources, which include star clusters, bright 
individual stars, and background galaxies, in the V 
band image using the IRAF task DAOFIND. We de¬ 
tect ^383,000 total sources and measure their bright¬ 
nesses within circular apertures varying from 0.5 to 5 
pixels with the background estimated within annuli of 

® Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science Insti¬ 
tute (STScI/NASA), the Space Telescope European Coordinating 
Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre 
(CADC/NRC/CSA). 
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Figure 1. Ground-based optical image of MlOl showing the lo¬ 
cation of all 10 HST ACS/WFC fields used in this work. The 5' 
scale bar is equivalent to ~9.3 kpc. 

7 to 13 pixels using the PHOT task within IRAF. We 
determine empirical aperture corrections out to 10 pixels 
from the curves of growth measured for 25 isolated clus¬ 
ters, and also apply an additional —0.107, —0.092, and 
—0.087 mag correct ion from 10 pixels to infinity for B, V, 
and I respectively (|Sirianni et al.l[2005l) . These aperture 
corrections are added to the measured photometry to ob¬ 
tain instrumental magnitudes. The apparent magnitude 
(mv) for each source in the VEGAMAG system is found 
by applying the following zero points: f43 5LF = 25.791 , 
F5551F = 25.738, and F8UW = 25.533 (IBohlinl 1200^ 
IMack et al.l[2nnl ]^ 

2.2. Cluster Selection 

In order to separate ancient star cluster candidates 
from our full catalog, we made a series of automated se¬ 
lective cuts followed by visual inspection of each object. 
We use the following automated selection criteria: 

1. Clusters brighter than my < 24.75 where my was 
measured within 3 pixel apertures (i.e. no aper¬ 
ture correction yet applied) to ensure high signal- 
to-noise ratio. 

2. Concentration index (Cl) > 1.15 where Cl is the 
difference between my measured within 1 pixel and 
3 pixel apertures to eliminate point sources; point 
sources have Cl values that peak around 1.00 with 
a standard deviation of 0.06. 

3. Colors, i.e. 0.55 < B - F < 2.0 and 0.75 < F - / < 
2.0, similar to those of Galactic GCs. 

Because Cl is a crude measure o f object size , we also 
use the BAOlab/ISHAPE software (|Larsenlll999f ) to mea¬ 
sure the EWHM of each object. ISHAPE fits proHles to 

^ http://www.stsci.edu/hst/acs/analysis/zeropoints 
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each candidate source to determine its FWHM (along 
with other parameters such as ellipticity). The light pro¬ 
file is a convolution of the point spread function (PSF) 
with a user determined function representing the spread 
in light from a cluster’s non-point-like size. We deter¬ 
mined the PSF of each of our fields by visually selecting 
^40 -50 isolated stars in each field. We choose a King pro¬ 
file (|K indll962|l with a concentration parameter (ratio of 
the tidal radius to the core-radius) of 30 to represent the 
cluster-like light profile convolved with the PSF for the 
ISFIAPE fitting. 

All objects with a measured FWFIM < 0.2 pixels 
(ce// « 0.46 pc) from ISFIAPE were removed from 
the catalog to further eliminate stars. We then plot¬ 
ted FWHM versus Cl and determined a three piece-wise 
linear fit to the relationship between Cl and FWHM. 
Upon visual inspection, objects falling outside of a per¬ 
pendicular distance of ^ 0.1 from the line only consisted 
of contaminants (usually faint, distorted “patches” that 
cannot be classified, noise, crowded point sources, etc). 
Therefore, only objects within a perpendicular distance 
of 0.1 from the CI-FWHM fits were kept. 

Finally, obvious background galaxies, chance superpo¬ 
sitions, and other contaminants were eliminated via vi¬ 
sual inspection. Figure [5] shows a few typical radial pro¬ 
files examined during the visual inspection of six selected 
candidate clusters and five stars (not included in the cat¬ 
alog). Note the clear difference between stars and the 
relatively well-resolved clusters. 

The final catalog consists of 326 candidate clusters (see 
Fig.[3]and Tabled]) with magnitude, color, and size mea¬ 
surements (six of these do not have FWHM measure¬ 
ments due to ISHAPE fit errors). For simplicity, we re¬ 
fer to the cluster candidates as “clusters” throughout the 
rest of the paper. 



Radius (pixels) 


Figure 2. Radial profile showing the enclosed percentage of fiux 
within a given radius (in pixels) for three faint clusters (dashed 
lines), three bright clusters (solid lines), and five stars (dotted 
lines). The point sources have significantly steeper profiles. 


2.3. Completeness 

To evaluate the completeness of our catalog of clusters, 
we generate artificial clusters, add them to the MlOl 
images, and re-run the detection and cluster selection 
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Figure 3. BVI color images of typical clusters from faint (top 
panel) and bright (bottom panel) cluster groups. Each panel is 
3.7”, or ~110 pc on a side. 


methods described in §2.1 and §2.2. We generate 4000 
artificial clusters using the BAOLAB task MKCMPPSF, 
which convolves a PSF with a user defined function, 
in this case KING30 profiles with two different input 
FWHM: 1.0 and 2.0 pixels (the motivation for these sizes 
is given in §3.3). We then use MKSYNTH to randomly 
place the artificial clusters in one of the images, where 
the magnitude range matches that of the real clusters. 
We detect sources with DAOFIND using the same pa¬ 
rameters as for the real clusters, measure photometry 
with PHOT, size measurements with ISHAPE, and run 
the automated cluster selection criteria and CTFWHM 
relation cut (see §2.2). 

Figure |4| shows the completeness fraction as a func¬ 
tion of apparent magnitude. Both sizes of artificial clus¬ 
ters show a decline in completeness as brightness de¬ 
creases with the more diffuse clusters declining at slightly 
brighter magnitudes. Approximately 80% of the artifi¬ 
cial clusters brighter than my = 23.0 make it through 
our selection pipeline, and approximately 50% of those 
brighter than my ~ 23.6. 

Although the clusters are randomly placed in the MlOl 
images, we ensured that 2000 are located within galacto- 
centric distance rgc ~ 2 kpc and 500 within rgc « 500 pc 
to ensure that the completeness of the inner regions of 
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Table 1 

MlOl Old Star Cluster Catalog^ 


ID 

02000 (hms) 

<52000 

V 

B-V 

V-I 

B - I 

reft (pc) 

945 

14 3 16.68 

54 19 38.04 

-5.89 

0.68 

1.08 

1.77 

5.3 

3072 

14 3 22.28 

54 19 50.40 

-5.76 

0.64 

0.97 

1.61 

3.8 

4335 

14 3 17.58 

54 19 55.68 

-7.35 

0.72 

0.98 

1.70 

3.6 

6166 

14 3 18.09 

54 20 3.09 

-7.05 

1.28 

1.75 

3.02 

2.4 

9429 

14 3 15.44 

54 20 16.08 

-6.77 

0.93 

1.34 

2.27 

2.0 

9936 

14 3 16.31 

54 20 18.48 

-5.96 

0.79 

1.16 

1.95 

7.3 

10101 

14 3 18.62 

54 20 19.33 

-6.30 

0.68 

1.04 

1.72 

7.6 

10379 

14 3 23.53 

54 20 20.59 

-6.49 

0.76 

1.21 

1.97 

4.7 

10633 

14 3 25.18 

54 20 21.79 

-8.83 

0.80 

1.19 

1.99 

2.3 

10654 

14 3 20.35 

54 20 21.89 

-5.89 

0.84 

1.19 

2.03 

7.7 


^ This table is available in its entirety in machine-readable form. 


MlOl are thoroughly tested. Figure [5] shows the com¬ 
pleteness fraction as a function of While the com¬ 
pleteness is slightly better for more compact clusters, 
there is no significant change in the completeness fraction 
with Tgc except within the innermost ^ 300 pc (dotted 
line in Fig. [5]). 


My 

-5 -6 -7 -8 -9 -10 



Figure 4. Fraction of selected artificial clusters versus ruv for 
FWHM= 1.0 (solid line) and 2.0 (dashed line) pixels. 

3. RESULTS AND ANALYSIS 
3.1. Cluster Luminosity Distribution 

Figure [5] shows the magnitude and luminosity distri¬ 
butions for the MlOl clust er sample alon gside those for 
MW GCs (taken from the iHarrla (|1996D catalog). The 
top panel shows a histogram with equal size magnitude 
bins, while the lower panel shows the luminosity distribu¬ 
tion with an equal number of clusters (seven) in each bin. 
In the magnitude distribution, the MlOl clusters have a 
similar shape to those in the MW (see §4 for further dis¬ 
cussion) at the bright end, but a drastically different one 
at the faint end. 

In the rest of this section, we compare the results of 
fitting different portions of the cluster luminosity distri¬ 
butions. These fits are performed with a simple least 
squares regression in logN-mv space. The slope, a, from 
these fits can be converted into the true power law index 
a for the LF dN{Ly)/dLy oc Ly by a = —(2.5a -I- 1). 
We first compare the results of fitting the bright end 
of the MlOl and MW GC distributions (see the bot- 



Figure 5. Fraction of selected artificial clusters versus rgc for 
FWHM= 1.0 (solid line) and 2.0 (dashed line) pixels. The dotted 
line represents the innermost region excluded from the fits in Figure 
1151 (see §3.4), which is the only significant drop in completeness as 
a function of distance. 

tom panel of Fig. HI). The fit range is My < —7.5, ap- 
proxim ately 0.2 mag bright er than the peak in the MW 
GCLF (| Ashman et al.l[T995r) . The best fits shown in Fig¬ 
ure H] are aMwi,bright = -1.97 ± 0.14 and auw,bright = 
— 1.91 ± 0.11. The standard deviation for 16 different 
binnings and ranges are ^ 0.16 for both the MlOl and 
MW data, comparable to the uncertainties. 

For the faint clusters, we use a magnitude range of 
—5.4 > My > —6.5. Here, completeness is a factor, and 
fainter than My ks —5.4, the LF declines steeply, indicat¬ 
ing the probable completeness limit, which matches well 
with the magnitude at which the completeness fraction 
of the artificial clusters drops below 50% (see Fig. 2]). 
We fit the MW GCs in the full faint range My > —6.54 
{my > 22.5 ) since the co mpleteness limit is not a con¬ 
cern for the iHarrisI (|1996r ) catalog. The slopes for these 
fits are crjv/ioi./ami — 2.03 i 0.05 and crjv/iv./ami — 

-0.46 ±0.170. The standard deviation for 16 different 
binnings and ranges (including more conservative com¬ 
pleteness limits brighter than My k, —5.4) are ^ 0.25 
and ^ 0.22 respectively. The shapes of the faint ends 
of the MlOl and MW GC distributions are very differ¬ 
ent. Note that any faint clusters not identified in our 
MlOl sample due to incompleteness will only steepen 

® Note that the value for a for the faint MW GCs is negative, 
while a, which is shown in the lower panel of Fig. [6] is positive. 
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Figure 6. Top: magnitude distribution (top) for our MlOl clus¬ 
ter sample (solid line), MW GCs (dotted line), and the difference 
between the MlOl clusters and the MW GCs (dashed line). A 
power law fit (dashed line) to the dashed line histogram for the 
bins brighter than the completeness limit gives a = —2.60 ib 0.26. 
The solid vertical line represents mv = 22.5 (My = —6.54), which 
is where we choose to divide our sample into “faint” and “bright” 
clusters (see text). The dotted vertical line represents the peak of 
the MW GC distribution. Bottom: luminosity distribution for 
our bright clusters (crosses), faint clusters (solid circles), and MW 
GGs (asterisks) with constant number binning. Fits to subranges 
of each group are shown with the same line styles as the top panel. 

OiMim faint and increase the discrepancy between the 
MlOl and MW GC distributions. 

Statistical tests confirm the fit results for the luminos¬ 
ity distribution of MlOl, that the bright end is similar to 
that in the MW, but the faint end is quite different. We 
compare the shapes of the MlOl and MW luminosity dis¬ 
tributions using the two-sided Kolmogorov-Smirnov (KS) 
test and the two-sided Cramer-von Mises (CvM) test. 
Both tests accept the null hypothesis (p-values > 0.1) 
that the MlOl luminosity distribution and the MW lu¬ 
minosity distribution are drawn from the same distribu¬ 
tion at the bright end My < —6.54, and very strongly 
reject the null hypothesis (p-values < 0.01) at the faint 
end My > —6.54. 

In fact, the fits do not produce p-values < 0.05 (strong 
rejection) until the samples include clusters fainter than 
My > —6.54. Therefore, we divide our sample into two 
groups throughout the rest of the paper, where “bright” 
clusters have My < —6.54 {my < 22.5) and “faint” clus¬ 
ters have My > —6.54 {my > 22.5). There are 90 clus¬ 


ters in the bright group and 236 clusters in the faint 
group. 

We perform a final fit to the faint clusters after sub¬ 
tracting the MW GC histogram from the MlOl cluster 
histogram (see the dashed line in the top panel of Fig¬ 
ure ED. No normalization of the MW GC histogram is 
performed prior to subtraction because the bright ends 
of the distributions are well matched (see Fig. ED ■ We find 
a = —2.60± 0.26 (uncertainties from Poisson errors, i.e. 
l/\/]V) in the range my = 23.58 to my = 22.03. Nine 
different combinations of bin sizes and data ranges (al¬ 
ways brighter than the completeness limit, my ~ 23.65) 
give a median a = —2.57 and a standard deviation of 
- 0.18. 

3.2. Cluster Colors and Luminosities 

In this section, we compare the colors and luminosities 
of old star clusters in MlOl with those of cluster popu¬ 
lations in other galaxies. The ages and metallicities of 
clusters older than ^ 1 Gyr become degenerate, making 
it difficult to establish differences in age based only on 
integrated broad-band colors; the colors are a stronger 
indicator of metallicity than of age in this regime. 

For bright clusters, the median B — V, V — I, and 
B — I colors are 0.72, 1.10, and 1.81, respectively, and 
for the faint clusters they are 0.65, 1.10, and 1.76. Fig¬ 
ure [3 shows B — V versus V — I colors with two metal¬ 
licity simple ste llar population (SSP) try ks (solar and 
Z = 0.008) from iBruzual fc Chariot! (120031 1 updated 2006 
data (from private communication). The cluster colors 
are approximately centered on the SSP tracks with some 
spread in the colors. This spread is almost certainly 
due to photometric errors, based on a comparison be¬ 
tween input and measured magnitudes for artificial clus¬ 
ters. Furthermore, the spread in cluster colors around 
the SSP tracks increases for fainter clusters, as expected 
when photometric errors are the dominant source of un¬ 
certainty. Therefore, the small differences between the 
median values for the faint versus bright cluster colors in 
B — V and B — I are unlikely to be strongly significant. 

The B — I color has the largest wavelength baseline 
and is therefore best suited to revealing multiple peaks 
in the colors, if they exist. Figure E] shows B — I his¬ 
tograms. The top plot shows the bright and faint MlOl 
cluster groups, and the bottom plot divides the clusters 
by distance, eith er inside or outsid e of 3 kpc from the 
center of MlOl ([Evans et al.ll201^ . The faint clusters 
drive this peaked distribution; the bright clusters appear 
to have a more even-color distribution with no obvious 
peak. Overlaid (dotted lines) are the mean blue and 
red peak B — I colors of eight brightest c luster galaxies 
(BCGs) determined bv iHarris et al.! (120061 1. Note that we 
have not converted from the VEGAMAG system to the 
Jo hnson-Gousins system . We estimate from Figure 21 
in iSirianni et al.l (120051) that with a conversion to the 
Johnson-Gousins system, the MlOl B — I colors would 
shift blueward by at most ^ —0.03. It is then clear that 
the total cluster histogram strongly peaks close to the 
blue BGG peak, with a weak tail of clusters extending to 
the red BGG peak. This red tail may be caused by red¬ 
dening within the host galaxy itself for which we have 
not accounted. In the bottom panel of Figure El the 
clusters within 3 kpc of the center of MlOl are slightly 
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Figure 7. Color index plot showing B — V vs. V — I colors for our 
full MIDI cluster sample faint (solid circles) and bright (crosses) 
groups and the MIDI YMCs with spectra (p luses). The lines show 
the evolution of IBruzual &; CharloU 1120031) (updated 2006 data) 
SSP models with solar (solid line) and Z = 0.008 (dashed line) 
metallicities. The 1 Gyr age on each SSP track is indicated by an 
open diamond. The MlOl clusters follow both SSP tracks well with 
some spread most likely due to photometric uncertainties. The 
arrow represents typical Galactic reddening for Ay = 0.5. Galactic 
reddening in the dire ction of MlOl is very low {E{B — V) < 0.01, 
IChandar et al.l 1120041 ^'). and we do not correct for it. 


redder than those outside of 3 kpc, although they do not 
have a peak centered on the red peak of the BCGs. The 
redder color could result from higher extinction or pos¬ 
sibly higher metallicity of these more centrally located 
clusters, but we cannot distinguish between these two 
possibilities with the currently available data. 

Figure [9] shows the CMD for our cluster sample 
with MW GCs, Large Magellanic Gloud (LMG) inter- 
medi ate age clu s ters (1-3 Gyr, magn itudes and colors 
from Bic a et al.l (119961) and a ges from iMucciarelli et ^ 
(12007); Piatti et al.l (|20^: Goudfrooii et al.l ( ^lll) : 
iPalma et al.l (|2013D . and iKerber et al.l (I2007D 1. and old 
open MW clusters (1 - 8 Gy r, magnitudes, colors, and 
ages from iLata et (j2002D l included for comparison. 
Again, note that a conversion to the Johnson-Gousins 
filters has not been applied to the MlOl clusters; we 
estimate B — V colors would shift redward by at most 
^ -1-0.04 (iSirianni et al.l[200^ . The MlOl cluster sample, 
as a whole, has colors similar to the MW GGs, but with 
significantly more clusters fainter than My = —6.54. It 
is of note that the MW GG sample plotted here is miss¬ 
ing B — V for 40 clusters, 34 of which are in the faint 
cluster region. Even taking this into account, however, 
there are still many more MlOl faint clusters (236) than 
MW GGs (53) in this region. 

The LMG intermediate age clusters also lie largely in 
the region of the GMD dimmer than My = —6.54, but 
within the B — V color range of the MlOl clusters. Old 
open clusters in the MW occupy a similar color space as 
the LMG clusters, but with even fainter magnitudes. The 
median colors for the faint MlOl clusters, intermediate 
age LMG clusters, the old open clusters, and MW GGs 
fainter than My = —6.54 are B — V = 0.65,0.67,0.71, 
and 0.74 respectively, with the MlOl clusters shifting up 
to -1-0.04 with conversion to the Johnson-Gousins system. 
All of these groups may have consistent B — V colors as 
the MW GG median value is affected by the 34 missing 
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Figure 8. Top: B — I color histogram with MlOl clusters divided 
into bright (dashed-dotted line) and faint (dashed line) groups. 
Vertical dotted li nes are the mean b lue and red peak B — I colors 
for eight BCGs JHarris et al.lf^06lh Bottom: same as above, 
except the clusters are now divided by distance with the dashed 
line representing clusters outside of 3 kpc from the center of the 
galaxy and the dashed-dotted line representing inside of 3 kpc. 
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Figure 9. B — V color magnitude diagram showing our MlOl 
sample (solid circles), MW GCs (asterisks), LMG intermediate age 
clusters (1-3 Gyr, open squares), and old open clusters (1-8 Gyr, 
open triangles). The dotted line corresponds to the luminosity 
division at My = —6.54 imposed on our MlOl sample (see §3.1). 
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values. 

Also shown in Figure 0 are the colors for 25 young 
massive clusters (YMCs, ~100s Myr) that are discussed 
further in §4.1; they are part of a spectroscopic study 
on clusters in MlOl (to be released in a follow-up pa¬ 
per). Their categorization as YMCs is derived from the 
strength of the Balmer lines seen in Gemmz/GMOS spec¬ 
tra. The YMC colors shown in Figure [3 however, are 
measured from the HST BVI images studied here with 
the same treatment as the rest of the cluster catalog. It 
is clear that the YMCs overall are much bluer than the 
faint and bright populations with a median B — V = 0.17 
and Y - / = 0.52. 

3.3. Sizes 

Figure [To] shows the histogram of cluster sizes (reff) for 
the bright and faint cluster subsamples. The FWHM of 
each cluster was measured using the ISHAPE software, 
as described in §2.2, then co nverted to r^ff assuming a 
distance to MlOl of 6.4 Mpc (iShapoee fc Stane^ 1201 111 . 
The figure shows that the median rgff values for the bright 
(2.41 pc) versus faint (4.27 pc) clusters are different, with 
fainter clusters tending to have larger sizes. MW GCs 
have a median reff of 2.98 pc, similar to that found for 
the bright clusters in MlOl. 

Figure |TT] plots cluster luminosity versus reff for the 
bright and faint MlOl samples, and for MW GCs. The 
MW GCs cover a similar parameter space as those in 
MlOl, with the bright clusters in MlOl being more com¬ 
pact with less scatter to larger radii than the faint clus¬ 
ters. 



Figure 10. reff or half-light radius histogram for our MlOl sample 
divided into bright (solid line) and faint (dashed line) clusters with 
the MW GCs (dotted line) also plotted. The median values of the 
distributions are shown by the vertical lines at the top of the plot 
with their respective line styles for each group. The median r^ff 
of the bright MlOl clusters (2.41 pc) is similar to that of the MW 
GCs (2.98 pc) while the faint MlOl clusters have a much larger 
median size (4.27 pc). 

In order to investigate whether or not there are system¬ 
atic biases in the sizes measured by ISHAPE for bright 
versus faint clusters, we compare ISHAPE measurements 
of artificial clusters of varying input sizes added to the 
MlOl images (see §2.3 for details). We use 2000 of 
the artificial clusters from the completeness testing with 
EWHM = 1.0 and 2.0 pixels in addition to 2000 artificial 



Figure 11. Absolute magnitude vs. reff for our MlOl catalog 
(solid circles) and MW GCs (asterisks). The horizontal dotted line 
represents the luminosity division at My = —6.54 imposed on our 
sample (see §3.1). The density of the faint MlOl clusters is much 
greater than that of the faint MW clusters. (Although the MW 
clusters are a more complete sample down to fainter magnitudes.) 
The brighter clusters have less spread in sizes and skew toward 
being more compact for both the MW and MlOl. 

clusters with EWHM =0.5 and 4.0 pixels. 

Figure [13 plots the cluster luminosity versus measured 
Teff for the MlOl clusters and the FWHM= 1.0 and 
2.0 pixel (Peff = 2.3 pc and 4.6 pc) artificial clusters. 
While it is apparent from the artificial clusters that there 
is a spread in the measured VeS which increases at fainter 
magnitudes, it is also apparent that this spread cannot 
fully describe the number of real, faint clusters with large 
sizes. 

We compare the results from artificial clusters with the 
observations using the Hodges-Lehmann (HL) estimator, 
which determines the shifts between the location param¬ 
eters of two data sets, in order to quantitatively establish 
whether the observed difference in cluster sizes between 
the bright and faint subsamples is real or the result of 
systematic and random errors. The results from the sta¬ 
tistical tests and artificial cluster experiments described 
below suggest that there is a physical difference in the 
sizes of bright and faint clusters in MIDI. 

The HL shift estimates the difference between the me¬ 
dian values of two data sets if the data sets are sym¬ 
metric about the median, or the difference between the 
pseudo-medians if the data sets are not symmetric about 
the median. It does this by computing the median of 
the differences between each pair of the values in the 
two data sets. It also gives confidence intervals and p- 
values for the significance of the separation as part of 
the Wilcoxon rank sum test. The HL estimator indi¬ 
cates a shift between the median sizes of bright versus 
faint clusters in MIDI of 1.26jlQ |^ ( errors are 95% confi¬ 
dence intervals). For the same magnitude ranges, there 
is a smaller shift of 0.30lj];);| (0.55t[];^f) for bright ver¬ 
sus faint artificial clusters with input FWHM values of 
1.0 pixels or reff = 2.3 pc (FWHM= 2.0, reS = 4.6 pc). 
We also tested the HL estimators for bright versus faint 
clusters for two faint binnings that excluded the faintest 
clusters and therefore highest uncertainty FWHM mea¬ 
surements (faint bins from —5.6 > My > —6.54 and 
—6.1 > My > —6.54; bright bins still My < —6.54). 
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Figure 12. The same as Fig. 1111 for the MlOl clusters (large 
solid circles) and artificially generated clusters (small solid circles). 
The vertical solid lines are the input sizes of the artificial clusters 
(top FWHM= 1.0 and bottom FWHM= 2.0). In the top plot, the 
real data are reasonably well matched by artificial clusters at bright 
magnitudes, but at faint magnitudes, it becomes apparent that too 
many real clusters have large sizes that cannot be explained by the 
spread due to measurement errors alone. In the bottom plot, we 
see that the faint clusters are better described by larger input size 
artificial clusters. 

This only increased the difference between the real and 
artificial cluster shift estimators. 

P-values for the real and artificial clusters for all bin- 
nings show that the shifts are significant with values 
ranging from 6.2 x 10“^^ to 0.03 (strongly accepting 
the alternative hypothesis that the true location shift 
is not equal to 0). We conclude that there is a physi¬ 
cal difference between the sizes of bright and faint clus¬ 
ters in MlOl, with the latter tending to be larger (by 
- 0.71 - 0.96 pc). 

3.4. Spatial Distribution 

The locations of our clusters within MlOl in R.A. and 
decl. are shown in Figure [T^ FigurelTTlshows the rgc dis¬ 
tribution for the faint versus bright clusters. The bright 
clusters are largely found within 9 kpc of the center, and 
quite centrally concentrated, while the faint clusters are 
more evenly distributed. 

Figure dni shows the surface density plots in logN-rgi‘^ 
and InN-rgc space. N for the surface density plots is 
derived by counting the number of clusters within the 



211.1 211.0 210.9 210.8 210.7 210.6 210.5 

RA(J2000) 

Figure 13. Positions in R.A. and decl. for our cluster sample. 
The bright clusters (crosses) are more centrally concentrated than 
the faint clusters (solid circles). 



Figure 14. Histogram of rgc for each cluster in our sample. The 
bright clusters (solid line) are concentrated toward the center while 
the faint clusters (dashed line) have a broader distribution. 

annulus inner and outer radii covered by the rgc bin. N is 
then divided by the area of the annulus to determine the 
cluster surface density. Here, rgc and N do not include 
the innermost 10" (^ 310 pc at the distance of MlOl) 
of the galaxy since this region is too bright to detect 
clusters down to the same magnitude limit as the rest 
of the galaxy. Beyond this region, the background level 
drops sharply and incompleteness does not vary strongly 
with rgc. This can be seen in Figure [5] which shows the 
completeness fraction as a function of rgc. 

We fit a de Vaucouleurs law of the form logN = 

bwauc + Cvaucfgc"*^ and an exponential of the form InN = 
bexp + CexpJ'gc to each of the three data sets which are 
reported in Table [5] (errors are Poisson errors based on 
the count in N). The de Vaucouleurs surface bright¬ 
ness profile is typically used to model spheroidal coin - 
ponents and elliptical galaxies (|de VaucouleursI 119481) , 
while exponential profiles pro vide a good des cription of 
disk components of galaxies (jPattersonlll94nl) . We also 
convert the slopes of the de Vaucouleurs fits into effec¬ 
tive radii {Re)^ the radius within which half of the total 
cluster population lies i n projection, with the formul a 
Re = (—3.33/cyauc)^ ()de Vaucouleurs fc Butal 119781) . 
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Table 2 

Coefficients for Surface Density Fits 


Data Set 

^Vauc 

^Vauc 

Re (kpc)'' 

^exp 

^exp 

Ro (kpc)*' 

Bright MlOl 

2.91 ±0.21 

-2.47 ±0.15 

3.28 ±0.80 

0.68 ±0.16 

-0.39 ±0.02 

2.57 ±0.13 

Faint MlOl 

2.17 ±0.19 

-1.70 ±0.11 

14.78 ±3.81 

0.64 ±0.14 

-0.25 ±0.02 

3.97 ±0.32 

MW GCs 

3.07 ±0.17 

-2.46 ±0.12 

3.35 ±0.66 

1.05 ±0.13 

-0.38 ±0.02 

2.61 ±0.14 


^ Effective radii calculated with the formula Re = (—3.33/cv^aiic) ^ llde Vaucouleurs fc Butairi978I L 
^ Scale radii calculated with the formula Rq = — l/cea:p. 


4.1. Are There Two Populations of Old Clusters in 
MlOl? 
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Figure 15. Surface density plots showing the area normalized 
number of clusters at different binned radii from the center of the 
host galaxy in log (left) and InN-r^c space (right). The 

symbols show the mid point of the radii bin they represent with 
either connecting lines to guide the eye (top) or best fit lines over 
laid (bottom). In all plots, the bright MlOl clusters (crosses with 
solid line) follow a centrally concentrated distribution, similar to 
the MW GCs (asterisks with dotted line), while the faint MlOl 
clusters (solid circles with dashed line) are much less centrally con¬ 
centrated. Slopes of the best fits to the data are shown in Table [2] 


The slopes of the exponential fits are converted into scale 
radii (i?o) by i?o = -l/cexp- 

Both the Re and Rq are almost identical for the bright 
MlOl clusters and MW GCs whereas the values for the 
faint MlOl clusters differ significantly from those of the 
MW GCs. Furthermore, the faint MlOl clusters do not 

appear to follow a straight line in logN-r\if^ space (see 
the left panels of Fig. [13, and so may not be well repre¬ 
sented by the de Vaucouleurs law. However, they follow 
a straighter line path on the InN-rgc space surface den¬ 
sity plot, indicating an exponential law better describes 
the distribution of the faint clusters. The bright clusters 

appear to follow a slightly straighter line in the logN-rgi'^ 
space than InN-r^e space, as expected from the high de¬ 
gree of steepness in linear space, and hence, the de Vau¬ 
couleurs law gives a better fit. 

4. DISCUSSION 

In this section, we first examine the properties (lumi¬ 
nosity, color, size, and spatial distributions) of the bright 
versus faint clusters to determine if our sample consists 
of two distinct populations of old star clusters. Then we 
more closely examine the properties of the faint clusters 
and compare to those in other galaxies. 


The shape of the LF gives the first indication that there 
may be two populations of red clusters in MlOl. As seen 
in Figure 13 the LF rises nearly continuously from the 
bright end down to the completeness limit, but shows a 
dip near My ^ —7.5 to —6.54. This dip can be explained 
if the LF is a combination of a peaked distribution (as 
seen for old GGs) at the bright end and a rising, power 
law-like distribution (as seen for younger clusters) at the 
faint end. For old GG systems in the MW and other 
galaxies, both luminosity and mass distributions have 
similar, peaked shapes, believed to result from mass-loss 
over ^ 12 Gyr, mostly driven by the effects of two-body 
relaxation. 

When we compare the MW GCLF to that of our MlOl 
cluster sample (see Fig. [Bland §3.1), we find the shapes 
quantitatively match well at the bright end, at least down 
to My = — 7.5, 0.2 mag brighte r than the peak in the 
MW GCLF (lAshman et al.lfl^ . The CMD (see Fig. [9]) 
shows that the colors of the majority of the bright clus¬ 
ters are similar to those of MW GCs. The median reff 
of the bright clusters is similar to the median size of the 
MW GCs (but smaller than that of the faint clusters; 
see §3.3). The bright MlOl clusters are centrally concen¬ 
trated toward the center of the galaxy, as are the MW 
GCs (see §3.4), which is expected for a spherically dis¬ 
tributed population. 

The similarity between the luminosity, magnitude, 
color, size, and spatial distributions of the bright MlOl 
clusters and the MW GCs strongly suggests that the for¬ 
mer are most likely a typical old (~ 12 Gyr) population 
of GCs in MlOl. Furthermore, our results suggest the 
MlOl bright clusters are largely metal-poor, consistent 
with the galaxy’s nearly bulge-less morphology. FigurejB] 
shows that the bright cluster colors skew toward the blue 
peak of BCGs ([Harris et al.ll200'^ . Corrections to the 
photometry that put our clusters on the same photo¬ 
metric system as the BCG data would push the colors 
further blueward as would any unaccounted for redden¬ 
ing corrections. Thus, the few clusters near the red peak 
of the BCGs, most likely do not account for a separate 
metal-rich population of clusters as seen in more bulge- 
dominated/spheriodal galaxies. 

Turning to the faint clusters, their LF can be described 
by a rising power law with an average slope of a = 
—2.6 ± 0.3 (see §3.1 and Fig. |B])- This is similar to what 
is found for y oung cluster populations in other galaxies. 
iLarseni (|2002ll found that young star cluster LFs in six 
spirals follow a power law of dN (Ly) /dLy oc LZ where 
a ranges from ^ —2.0 to —2.6. [Whitmore et al'l (|2014[ 1 
examined the LFs of clusters in 20 nearby (4 — 30 Mpc) 
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star-forming galaxies, and found an average a = —2.37 
with rms scatter of 0.18. 

In checking for similarities between the faint clusters 
and the MW GCs, we find the median reff of the MW 
GCs is smaller than that of the faint clusters, and their 
spatial distributions are very different. While the MW 
GCs and bright red MlOl clusters are concentrated to the 
center of their host galaxies, the faint clusters in MlOl 
are more uniformly distributed as expected if they reside 
in the disk. We conclude that the faint, red clus¬ 
ters form a population distinct from the typical 
old GCs, one which is associated with the disk of 
MlOl. 

4.2. What Are the Faint, Red Clusters in MlOl? 

Discovering the true nature of the faint clusters de¬ 
pends on determining their ages, metallicities, and 
masses. It is clear from Figure [7] that exact ages and 
metallicities cannot be determined accurately for each 
cluster from optical colors alone since the SSP tracks lie 
on top of each other in this color space. Nonetheless, we 
can compare the colors, luminosities, and sizes of these 
clusters with those in other galaxies where the cluster 
properties are better known. 

The faint clusters may be a population similar to the 
LMC intermediate age clusters or the brightest old open 
clusters since they share a similar color-magnitude space 
(see Fig. [9]), although the faint MlOl clusters do ex¬ 
tend to brighter magnitudes th an the old open c l usters . 
The faint fuzzies disc o vered by iBrodie fc LarsenI (|20n2[ l 
and iLarsen fc Brodid (|200Clll in lenticular galaxies have 
V — I ^ 1.0 — 1.5 which is similar to that of the faint 
clusters, ^ 1.1. Despite the similar colors and magni¬ 
tudes between faint fuzzies and our faint clusters, only 
48 of the 230 faint clusters with FWHM measurements 
have reff > 7 pc, but the faint clusters do still have a 
large median size, greater than the bright clusters or MW 
GCs. 

We can use the observed power law shape of the LF 
to constrain the ages of the faint clusters. As mentioned 
previously, the peaked shape found for old GCs is driven 
by stellar mas s loss d ue to internal relaxation. Figure 3 in 
iFall fc Zhand (|2001h . shows how a peak develops at the 
low end of the mass function and then moves to higher 
masses as the population ages, while the high mass end 
continues to have a power law shape. Since we do not 
observe a peak in the LF of the faint, red clusters, it must 
occur below our completeness limit. When luminosities 
are converted to masses, the age at which the peak of the 
mass function is just below the completeness limit is the 
maximum age of the faint clusters. 

We make the following assumptions for our calcu¬ 
lations: a stellar-mass loss rate for each cluster of 
Hev = 10~^ M(T,vr~\ a typical value for MW GCs 
(e.g., iFall fc Zhand (|2001ll l. and a single age r for the 
cluster population. This gives a predicted peak mass, 
Afp = /TevT, at different ages. We then use the SSP 
model predicted, age-dependent M/Ly to convert our 
observed cluster luminosities to masses at different as¬ 
sumed ages, and determine the mass equivalent of our 
completeness limit, Mum. Based on this methodology, 
the maximum age for our clusters is constrained by the 
value where Mp = M\{^. 

We find the maximum age to be ^ 700 Myr 1 Gyr) 


for Z = 0.008 {Z = 0.02). Since the faint clusters 
are found in the disk, we might expect the maximum 
age based on the more metal-rich model to give a bet¬ 
ter estimate, but MW old open clusters span met allic¬ 
ities from just above solar down to [Fe/H]^ —1 (jFriell 
Il995fl. Likewise, faint f uzzies have typical [Fe/H]^ —0.6 
(IBrodie fc Larsenl2002il , and LM C intermediate age clus¬ 
ters have [Fe/H]^ —0.4 to —0.7 (iPalma et al.ll20i3ll . 

The faint MlOl clusters do have a lower density 
than typical MW GCs, and the mass-loss rate due to 
relaxation-driven evaporation depends on the internal 
density of the clusters, Ue-v oi pj , where Ph is the half- 
light density (|Chandar et al.l 120011 : iMcLaughlin fc FaUl 
l2008fl . Therefore, a better estimate of the maximum age 
of the faint MlOl clusters might be determined by scaling 
the mass loss rate to the median size of the faint clusters 
(Mev.MlOl = liev,Mw(Ph,M10l/Ph,Mw)^^^)- Now, we find 
the maximum age to be ^ 9 — 10 Gyr for Z — 0.008 and 
^12—13 Gyr for Z = 0.02. Thus it is possible that the 
faint clusters might be quite old. 

We can also put constraints on the minimum age of the 
faint clusters. iBarmbv et al.l (|2006[1 suggested that the 
faint, red cluster population in MlOl may be reddened, 
young disk clusters; however, they do not provide any 
estimated ages. We examine B — I color images which 
highlight the locations of dust lanes in MlOl and find 
no preference for the faint clusters to be embedded in 
dust. Likewise, we examine archival HST/WFPC2 Ha 
(F656N) images and find no preference for the faint clus¬ 
ters to be in or near H II regions. This strongly indicates 
that they are older than 10 Myr, as clusters younger than 
this are expected to have not yet fully dispersed, leaving 
the gas and dust clouds from which they formed. 

Also, we have Gemini/GMOS spectra for 25 YMC 
clusters which exhibit strong Balmer lines with no Ha 
emission lines, indicating ages of a few hundred Myr. 
We find that their median colors from our FIST data 
are B — V ^0.2 and H ~ 0.5 (see Fig. [T]), sig¬ 
nificantly bluer than the faint red clusters, which h ave 
B — V ^ 0.7 and V — I ^ 1.1. IBarmbv et al.l (|2n06[l se¬ 
lected 1260 “blue” clusters in MlOl with {B—V)o < 0.45 
and H < 23, similar to our spectroscopically confirmed 
YMCs, and they found these clusters appear to coincide 
with the spiral structure of the MlOl disk. The faint red 
clusters likely do not follow the spiral structure (despite 
being in the disk), because they have had sufficient time 
to disperse throughout the disk and away from their birth 
sites. We therefore believe that the faint red clusters are 
older than a few hundred Myr. 

The faint red clusters are associated with the disk 
of MlOl, just as old open clu sters are in the MW 
disk ljPortegies Zwart et al.l 12010( 1. Note that the scale 
radius from the exponential fit to the faint MlOl clus¬ 
ters (^ 4 kpc) is similar to the scale radius of the 
disk of MlOl according t o optic al photometry (i?o =4.6- 
4.8 kpc lOkamura et al.l (I1976IH . Interestingly, faint 
fuzzies have also been found to be as sociated with the 
disks of their lenticular host g a laxies (I Brodie fc LarsenI 
120021 : iGhies-Santos et aP 120131: llbrbes et aT 201411 de - 


spite their ages >7 — 8 Gvr. iChies-Santos et al.l (|2013ri 
find that the strong association of the NGC 1023 faint 
fuzzies with their galaxy’s disk-rather than its bulge-to 
be evidence that the faint fuzzies are simply very old 
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open clusters. They further predict that such clusters 
should be found in spiral galaxy disks as well, and that 
the only reason they have thus far only been identified 
in lenticular galaxies is because of their smooth disks, 
which make them easier to observe. The excellent reso¬ 
lution of HST /ACS allo ws us to identify such clust ers in 
MlOl for the first time. iScheepmaker et ^ (|2007f) also 
found six clusters in the spiral galaxy M51 that match 
the definition of faint fuzzies. We conclude that the 
faint clusters studied here are old (r > 700 Myr) 
and part of the disk of MlOl, similar to but more 
massive than old open clusters in the MW. 

5. CONCLUSIONS 

MlOl appears to have two populations of old star clus¬ 
ters: a typical population of old GCs and a fainter pop¬ 
ulation of intermediate-age to old disk clusters. For the 
population of old GCs, we find: 

1. Their luminosity distribution, colors, sizes, and 
spatial distribution are similar to those of the MW 
GCs. 

2. Their spatial distribution shows a central concen¬ 
tration which is consistent with a spherically dis¬ 
tributed halo population. 

3. Their colors are skewed toward the typical blue 
peak of BCGs which indicates they are most likely a 
metal-poor dominated population, fitting with the 
nearly bulge-less morphology of MlOl. 

For the fainter population, we find: 

1. Their luminosity distribution is similar to the 
power law shape of young cluster populations in 
other spiral galaxies, dN{Lv)/dLv oc Ly where 
typically a ~ —2 to —2.6. We find the average 
power law fit to our faint clusters in MlOl has a 
slope of a = —2.6 ± 0.3. 

2. Age constraints determined from the shape of the 
LF indicate that they could be quite old, up to 
~ 12—13 Gyr, although this constraint depends 
on a number of assumptions. 

3. They are older than a few hundred Myr because 
their colors are much redder than a sample of MlOl 
YMCs with spectroscopic ages of a few hundred 
Myr. These clusters also do not follow the spiral 
arm st ructure that the blue MlOl clusters identi¬ 
fied bv iBarmbv et al.l (|2006ll do. 

4. They have a fairly extended spatial distribution, 
quite different from the centrally concentrated, 
bright GCs. They are most likely associated with 
the disk. 

5. They occupy the same luminosity-color space as 
LMC intermediate age clusters, the brightest old 
open MW clusters, and faint fuzzies. Old open 
clusters and faint fuzzi es are also located in the 
disks of their galaxies. iGhies-Santos et aP (l2013f ) 
concluded that faint fuzzies are analogous to old 
open clusters and should be found in spiral galax¬ 
ies. We conclude the MlOl faint clusters are most 
likely these old disk clusters. 


The evidence for a large population of old disk clus¬ 
ters in MlOl shows that the peak of the GCLF may not 
be an accurate distance indicator for all galaxy types, 
especially not spiral and lenticular galaxies. We would 
caution against using it as a stand alone measure of dis¬ 
tance. Other spiral galaxies need to be examined for 
faint, disk cluster populations. 

We are in the process of using Gemini GMOS spec¬ 
tra of bright, old GCs in MlOl to measure metallicities 
for calibrating the color-metallcity relationship for MlOl. 
Such a relationship can be applied to our bright cluster 
sample to determine more accurate ages and metallic¬ 
ities than can be found from the optical colors alone. 
Unfortunately, we do not have spectra for any clusters in 
the faint population. Preliminary results show that the 
bright cluster subset are quite old, with ages similar to 
those of Galactic GCs. Results from this analysis will be 
presented in a follow-up paper. 
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